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We present a method to calculate i/i/ energy losses from neutron star crusts, which automatically 
takes into account for collective effects, and allows to calculate the total emissivity without a separate 
consideration of particular processes. We show that the formula we obtain describe the known results 
for the emissivity due to plasmon decay and Bremsstrahlung from degenerate electrons, when one 
of this processes dominates. In the case of low temperatures, our formula gives a suppression of 
the electron vector weak-current contribution to i^i? Bremsstrahlung, due to the collective effects 
. discussed in this paper. 
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I. INTRODUCTION 



According to modern theoretical scenarios, within the first hundred years of existence of a neutron star the cooling of its 
interior layers occurs in an extremely nonuniform way. Due to the so-called "fast cooling processes" in nuclear matter, 
the core of a neutron star cools down very fast, while the external crust cools down much slowly, by heat diffusion 
towards the inside and by neutrino radiation ^ . The time necessary for thermal relaxation crucially depends on the 
intensity of neutrino and antineutrino radiation from the crust of the neutron star, which influences the temperature 
of the stellar surface, and thus the gamma and X-ray which are observed. It is well-known that radiation of neutrinos 
[ and antineutrinos from the neutron star crust is caused dominantly by Bremsstrahlung of degenerate electrons, 
• and by plasmon decay into neutrino pairs. The Bremsstrahlung process, operating in a liquid or crystalline phase of 
I the crust, has been previously studied by many authors |^ - g. These calculations have been performed neglecting 
f~| . collective interactions of electrons with the neutrino field. When the temperature of the crust is smaller than the 
Oh' electron plasma frequency, such an approximation is not justified because, for the scenario under consideration, the 
Q , wavelength A of radiated neutrinos and antineutrinos is larger than the electron Debye screening distance Dg. By 
' undergoing a quantum transition, the radiating electron electromagnetically induces a motion of other electrons inside 
[ the Debye sphere around itself. The weak current of perturbed electrons inside the Debye sphere generates neutrinos 
d • coherently with the initial electron, therefore screening its vector weak coupling with the neutrino field. Since this 
^ [ collective effect takes place when the neutrino wave-length is of the order, or larger, than the Debye screening distance, 
•'H . the relevant parameter of the problem is k^D^, where k is the momentum carried out by the neutrino-paii^ The 
' screening effect was demonstrated in |^ for neutrino-pair emission due to electron-phonon scattering in a crystalline 
crust, where the condition fc^-Dg ^ 1 is fulfilled. In this limiting case, valid when the temperature T is much smaller 
than the electron plasma frequency LOpe, the effective vector weak-current of electrons is totally screened by the plasma 
polarization, which dramatically modifies the neutrino emissivity. 

In the present paper, we introduce a method of calculation for vi? energy losses from a neutron star crust, which 
incorporates collective effects in a degenerate plasma at arbitrary temperatures, limited only by the degeneration 
condition T <C /^e , where /ig is the electron chemical potential. This condition holds during the cooling epoch 
described above. 

Collective effects appear if one takes into account the possibility of virtual photon-exchange among electrons in 
the plasma. To obtain the neutrino emissivity, we use the fluctuation-dissipation theorem in order to relate the 
weak current-current correlation function to the imaginary part of the exact retarded polarization functions of the 
plasma. This procedure, as we will see, has the advantage that collective effects are automatically included, and 
do not need a separate consideration. The obtained formula for the neutrino energy losses includes contribution of 



^In what follows we use the system of units h = c — 1 and the Boltzmann constant ~ 1- The fine-structure constant is 
a = e^ = 1/137. 
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plasmon decay into neutrino pairs, as well as the process of neutrino-pair emission due to electron collisions with 
nuclei and phonons. We identify the latter mechanism as the Bremsstrahlung process from electrons. Thus, the 
total neutrino-pair emissivity of the crust, caused by vv decay of plasmons and Bremsstrahlung from electrons is 
represented by a unique expression. 

This paper is organized as follows. In Sect. 2 we discuss the reference status of the problem, and derive some 
formulae describing the rate of vv Bremsstrahlung, which we use in the following considerations. In Sect. 3, we show 
how the same vv Bremsstrahlung rate can be obtained by the use of the Optical Theorem, and calculate the imaginary 
parts of the retarded polarization tensors of a liquid or crystallin crust. In order to incorporate collective effects, in 
Sec. 4 we generalize the results obtained by the Optical Theorem by replacing the imaginary parts of the polarization 
tensors by the weak current-current correlation function, according to the Fluctuation-Dissipation Theorem. This 
correlation function is obtained by summation of all diagrams with allowance of intermediate photon-exchange among 
the electrons. In Sect. 5 we derive a general formula for the vv emissivity, which includes the energy losses due to 
electron Bremsstrahlung and plasmon decay. Some limiting cases and numerical tests are presented in Sect. 6 in order 
to demonstrate the validity of the obtained general formula. Discussion of the results and conclusions are shown in 
Sect. 7. 



II. STATEMENT OF THE PROBLEM 

We use the Standard Model of weak interactions, and consider low-energy electrons, which are typical for neutron 
star interiors. Therefore the in- vacuum weak interaction of electrons with the neutrino field can be written in a 
point-like current-current approach 

neff = ^J r{x)J,(x)d'x, (1) 

where Gp '^s the Fermi coupling constant, and 

= Prjt" (1 - 75) 1/ 

is the neutrino current. The vacuum weak current of an electron is of the standard form, which is a sum of vector 
and axial-vector pieces 

= ?/i (Cy7^ - Ca7p75) V'- (2) 

Here, tp stands for electron field; Cy ==5 + 2 sin^ 6w , Ca — \ stand for emission of electron neutrinos, whereas 
C'y = — ^ + 2 sin^ 9w , Ca — ~h ^^"^ ^'^ used for muon and tau neutrinos; 6w is the Weinberg angle. 
The differential rate of vv production takes the form 

2 (2^)^ J 2E J 2E'-' ^ ^ ^ 

Here AI^ is the matrix element of the weak transition current for the electron in the Bremsstrahlung process. The 
symbol Tr (M^Mj^) includes statistical averaging and summation over initial and final states of the background (see 
below). Summation over all initial p = {E, p) and final p' = {E' , p') states of the electron has to take into account the 
Pauli principle via the appropriate blocking factors, with the Fermi distribution function / (E) of degenerate electrons. 
Finally, ki — (cji,ki) , k2 — (102,^^2) are the neutrino and antineutrino four- momenta, respectively. To simplify the 
calculations in what follows, we consider an ultrarelativistic electron gas, which is typical for the dominant volume of 
the neutron star crust. 

The matrix element for the weak transition current of the electron has the following form 

Af^ = u (p') [7^ {Cv - CaIz) G {p' + X) 7" + 7"G (p - K) 7^ {Cy ~ CaI^)] u (p) 

/OC 
dte-'""' (/I A„ (-q, t) 6t (p', a') b {p, a) (fci) (fca) \i) . (4) 
-00 

In this equation, p — {E, p) and p' — (£", p') are the initial and final momenta of the electron. The operators b, c, 
d refer to electrons, neutrinos and antineutrinos, respectively. In this process, the momentum and energy transferred 
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to the background are q = p — p' — k and qo = E — E' — uj, where K = (w, k) is the total four- momentum of the 
neutrino pair : K = ki + k2. . Within the Coulomb gauge, used in this section, the electromagnetic potential has 
only a non- vanishing (scalar) component. 

The energy transfer to the background, qo, is small with respect to the momentum transfer |q|, and is therefore 
neglected in the electron matrix element. On the other hand, the momentum transfer to the background, of the order 
|q| ~ pf , where pp is the electron Fermi momentum, is large with respect to the neutrino-pair momentum, therefore 
we assume that q=p — p' — kc^p — p'. Keeping this in mind, we can use the well-known method of soft photons 
|13| . With the aid of commutation rules and the Dirac equation, one can represent the matrix element as follows 



u{p')j^ {Cv-CA75)uip) 



nOO 

X / dte-'** (/I Ao (-q, t) fot (p', a') b {p, a) ct (fci) dt (^2) |z) 



(5) 



After squaring and taking the trace, and performing summation over final states and thermal averaging over initial 
states of the background we obtain 

Tr {M^f^W) = 4 (C^ + Cl) j dqodPqS (p - p' - q - k) (5 (^ - - go - w) 



{2EE' - pp') 



p- 



{AoAo) 



(6) 



^p'K pK J \p'K pK 

The Fourier transform of the A A correlation function can be written with the aid of the dynamic form factor S (qo, q) 



(^oAo),,, = f^— (2^(90, q). 

(r£(q)) 

The dielectric function of a degenerate ultrarelativistic electron gas can be taken in the static limit 



e (q) = 1 + 
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q2i?2 ^3 



2 l-3x^ 



6x 
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where 



(7) 



(8) 



2pf 



Using the approximations discussed above, we can rewrite Eq. (^ as follows 

S{qo,<l) {2EE'~pp') 



Tr (Mt^Af'') =4{C^ + Cl) / dq^d^^q 



Vg2£ {q) 



(9) 



(10) 



In a liquid crust, the energy transfer to the background is negligible due to the large nucleus mass Mi. In this case we 
can neglect go everywhere except in the dynamic form factor. This allows the integration over dq^ to be performed. 
We obtain 



Tr (Mt^Af-),,^^.^ = STTiV, {Cl + C\) J d\ 



p'^ p^ 
yP'K pK 

The static structure factor of ions is defined as 



J \p'K pK J ^ 



Sst (q) {2EE' - pp') 
p-p' -q)S{E- E' -uj). 



(11) 



(12) 
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where Ni is the number density of ions. We use the structure factor calculated in [|5[ for a one-component classical 
plasma, which is a tabulated function of the non-ideality parameter F of the ionic component of the plasma. 
In the case of a crystallin crust, the dynamic form factor can be written, in the one-phonon approximation |H] : 



S (go,q) = 27r^e" 



-2W{q) 



E 

sAK 



(q ■ ^As) 
2wAs 



(5 (<70 - t^As) (^q.s-K , S {qo + UJXs) Sq-s-K 
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(13) 



Here ujxs is the phonon frequency, which depends on the wave vector s. Each phonon mode A is defined by its 
polarization vector e\s. The K 's are reciprocal lattice vectors. The Debye- Waller factor is of the form 



2W{q) 



(q ■ ^As) 

sAK 



coth 



/^As\ 

\2t) 



(14) 



III. TREATMENT BY THE OPTICAL THEOREM 

With the help of the Optical Theorem, the rate Eq. (^ of vv production can be written as 

dT=^ t Im [n,. {K) Tr {ff*)] , (o\, ' (15) 

2 exp (y) - 1 2u;i(2T:y 2u;2(2T:y 

The exact, irreducible retarded polarization tensor H^'^ of the medium, represents the sum of compact diagrams 
which include inside the electromagnetic interactions of the electron with nuclei, and have ends at the weak vertex 
(Cy7/i — C'aTmTs)- This tensor is the following sum of vector and axial pieces 

- ^""v + S^n-. (16) 

We traditionally include an extra-factor Ane^ in the definition of all polarization tensors. By this reason, the factor 
l/47re^ has been included in Eq. ([T^). Following these notations, the vector- vector tensor 11'^ is the retarded 
tensor for the electromagnetic polarization of the plasma. In the absence of external magnetic fields, the parity- 
violating axial- vector polarization does not contribute to the rate of neutrino-pair production. In fact, by inserting 
J d^KS^^^ {K — fci — = 1 in this equation, and making use of the Lenard's integral 

4:7T 

= — {K^K, - K^g^,) e (K^) {K°) , (17) 
where 9{x) is the Heaviside step function, Eq. (|l^) can be written as 



3 2 ..)_Ti-n-(i.)(K,i..-KV)^m^m(^. (18) 



dT 

Since the axial-vector polarization has to be an antisymmetric tensor, its contraction in (|18|) with the symmetric 
tensor K^K^ — K'^g^i, vanishes. 

A. Polarization tensors 

To specify the components of the polarization tensors, we select a basis constructed from the following orthogonal 
four-vectors 

h^^^^, (19) 

where the space-like unit vector n — k/fc is directed along the electromagnetic wave-vector k. Thus, the longitudinal 
basis tensor can be chosen as Lp^ = —IpW- , with normalization = 1. The transverse (with respect to k) components 
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of n'''^ have a tensor structure proportional to the tensor T^f^ = {gP^ — hPh^^ + Fl^^), where g^^ = diag(l, — 1, — 1, — 1) 
is the signature tensor. This choice of T'''* allows us to describe the two remaining directions orthogonal to h and 
I. Therefore, the transverse basis tensor has normalization TP — 2. One can also check the following orthogonality 
relations: IpT^^ = 0, as well as kiT'^^' = 0, and KpL^P — KpT^P = 0. In this basis, the vector-vector polarization 
tensor has the following form 



[K) = n (K) LPP + 7Tt [K) TPP, 



(20) 



where the longitudinal polarization function is defined as tti {K) = (l — cj^/fc^) Hy" and the transverse polarization 



function is found to be TTf [K) — {gppYiy 
They can be written as 



) /2. The axial-vector polarizations have to be antisymmetric tensor^. 



(21) 



where ePP^^ is the completely antisymmetric tensor (^^oi^s _ -kav {K) and i^vA (K) are the axial-vector polar- 
ization functions of the medium. As for the axial term, it must be a symmetrical tensor. The most general expression 
for this tensor is, therefore 



(K) = TTi (K) LP"" + TTt (K) TP" + TTA (K) gP". 



(22) 



Real parts of the retarded polarization tensors are the same as the real parts of time-ordered (causal) polarizations, 
which can be written in the one-loop approximation as 



n^'' = 47rie^Tr 



Y G{p) Y G{p + K) 



(23) 



li^P^ = Ame^ Tr 



d'^p 



YG{p)Yl^G{p + K) 



(24) 



U'Y = Ame^ Tr 



d^p 
(2^ 



ri^G{p)Yl5G{p + K) 



(25) 



Here, G{p) is the in-medium electron propagator, which includes the Pauli principle restrictions. This approximation 
has been studied by different authors. It corresponds to a coUisionless plasma, and thus the one-loop polarization 
functions are real-valued in the case > 0, when Landau damping is not possible. In an ultrarelativistic, strongly- 
degenerate electron plasma, the longitudinal and transverse polarization functions take the formQ [ pO| : 



1 / cj' 



TTT In — 

2k LU 



(26) 



Re TTt = -Wpg 



1-^ 



OJ 
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UJ , ^ - 
■ In ■ 



J \ 2kvp UJ — k 
The electron plasma frequency and the Debye screening distance are defined as 

9 AnUee'^ 1 



fJ-e 



JJ2 ~ ^^Pe' 



(27) 



(28) 



with — Pf and n^. being the chemical potential (the Fermi energy) and the number density, respectively, of 
degenerate electrons. The axial-vector one-loop polarization tensor is given by 



■^We consider also the axial- vector polarization because it will be used in the next Section. 

^Our Eq. (H) differs from Eq. (A39) of the Ref. @ by an extra factor {cj^/k^ - l) because our basis Z^, hP is different from 
that used by Braaten and Segel. All components of the complete tensor Eq. (^^ identically coincide with that obtained in jl^] 
for the degenerate case. By the same reason, an extra factor V appears in the tvva expression (H). 
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Re -KvA (K) = —pfViO f 1 - 



LU UJ + k 

1 — — In 

2k UJ — k 



(29) 



Imaginary parts of polarization tensors appear at next order of corrections, when electron collisions with ambient 
particles arc included. The partial contribution of the electron scattering off nuclei to the imaginary part of the 
retarded polarizations can be obtained by comparing the differential rate of the vV production Eq. ( [l^ ) with Eq. (|^). 
We arrive to 



exp(f) -1 



Im 



^2 



Tr [M^^'M") , 



(30) 



where summation over initial and final states of the electron should be understood as described above. This formula 
can be identified with the optical theorem. By using the basis expansion Eqs. (pO|-22) and the following relations 



(31) 



we obtain 



(32) 



(33) 



luiTTA = 2^e2!^4^-^ ^ h^h, Tr {M^^M'') 



^A 



(34) 



A direct evaluation of formulae Eq. d32-R3) yield, for the liquid crust 



Imni {uj,k) = ujl^Za^^Fi {u,T) 

TT K 



(35) 



ImTTt (w,fc) = --Ldl^Za^^^Fo (m, T) - -^ImTTj, 
■K ^ k 2 



Im-KA = 0. 

The functions Fi {u, T) and Fq {u, T) are defined by the following integrals 



(36) 
(37) 



f;'^^'' (u,r) 



1 



2pf 



dq Sst (g) 
\e{q)\\ 



1 



1 ^ m 

2 ux(, Q — xu 



(38) 



with X = q/2pp^ and C, = ^yl — u^Y^-x^. Since we consider point-like nuclei, the fimctions Fi {u,T) and Fq {u,T) 
depend only on the variable u = k/u and the non- ideality parameter of the ionic component of the plasma |ll|: 

72„2 75/3 „^ 

= 0.02254^^, (40) 
ai la nif. 



where a is the ion-sphere radius. Expressions Eq. (|3^) and Eq. ( |39D are valid for T < 172. Above this value, the 
crust crystallizes. For the crystallin crust we obtain the more complicate expressions : 
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1 



2pF 



dq I — 
NiUUJ J„ \e{q)fq x'^ 

exp (^) - 1 



1 — - ^ In 

2 uxQ Q — xu 



, (w + go) c-. X 



(41) 



1 



dq 1 — x"^ C + xu 



In- 



^jw7o \e{q)\^q "' ( - xu 

exp (^) - 1 



, + go) c, X 



exp 



1 



(42) 



Inserting the explicit form Eq. (13) of the one-phonon dynamic form factor, and performing the integration over go 
we obtain 



sAK 



2pF 



dq 6-21^(9) (q.e^^)l- x"" 1 - u\ C + xu 



1 



In- 



|£(g)l^9 V 2uxC,"'C-xu^ 

^(q-s + K) {uj + ujxs) , (5(q + s + K) {uj - ujxs) 



1 



— 1 exp ( 



- 1 



(43) 



sAK 



2pF 



dq e-^W'C?) (q . J 1 - x'^ ( + xu 



k(g)rg c-a:w 

(5(q- s + K) {uj + ujxs) , 5(q + s + K) (w - was) 
1-e 



(44) 



IV. INCLUDING COLLECTIVE EFFECTS 



By the Fluctuation-Dissipation Theorem, the function 

^^-^ (X) = Im n^'' {K) (45) 

exp (j.) - 1 

should be identified with the Fourier transform of the correlation function of two distant weak-currents in the plasma. 
However, the above approximation, with H^"^ being the irreducible polarization tensor, does not take into account 
for collective effects in the correlation function. To generalize the correlation function to this case, we should also 
include the exchange of an intermediate virtual photon between electrons. This can be done by summation of the two 
diagrams shown in Fig. 1, where the thick dashed line is the in-medium photon propagator D^^^ {K), defined as the 
infinite sum of graphs shown in Fig. 2. In this figure, the thin dashed- line represents D^^ (K) - the photon propagator 
in vacuum. According to this series, the exact photon propagator in the medium satisfies the Dyson's equation 

D,, (K) = D% {K) + j-Dl^ [K) n7 {K) D^, {K) (46) 

Within the Lorentz gauge, Tf^v i-^) ^^e following form 



4tt 

D°Ap (K) = {gxp - KK) , (47) 
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The compact block Ily' {K) is the retarded polarization tensor of the plasma, given by Eq. (|2C|). The retarded 
propagator of the in-medium photon has the same tensor structure as the vector-vector polarization tensor. The 
solution to the Dyson's equation is of the form: 



Dxp [K] = Di (K) Lxp + Dt (K) Txp 



(48) 



with 



A {K) = 



47r 



- TTl (K) 

Thus, the sum of graphs of the Fig. 1 gives the new tensors : 

1 



'In 



~ TTt {K) ■ 



ii>^' (k, oj) - u'^^ (k, Lo) + — n^.^ (k, u) Dxp (k, c.) n^r (k, c.) , 



1 



(k, = n>x^ (k, u) + (k, D^p (k, u) n^y\ (k, u) 



By substituting Eqs. <^^^ , (Hsj), and into Eqs. @ we obtain 



(49) 



(50) 



'^W — lp2 ^~ 1P2 ' 

ii^ — TTl A ^ — TTi 



(51) 



-TT^" 



(52) 



— TTt 



(53) 



Replacement of the irreducible tensors in Eq. (g^ by the generalized expressions ( pl[f33| ) results in the following weak 
current-current correlation tensor 



Im 



exp I 



47re2 \K^- 



T{2 

J^^^ _|_ ^_rppi 



CI 



47re^ V ^ — TTt 



(54) 



V. ENERGY LOSS DUE TO NEUTRINO-PAIR RADIATION 

We consider the emissivity of neutrino pairs from the plasma , i.e. the total energy which is emitted into neutrino 
pairs per unit volume and time. Taking into account collective effects, to the lowest order in Gf, the emissivity is 
given by the following formula: 



Q 



2 3 {2nf^Jo 



duo uj \ d^k {K^K, - K^gp,) 1'^'' {K) . 



(55) 



The symbol 'Y^^ indicates that summation over the three neutrino types has to be performed, with the corresponding 
values of Cy and Ca, as explained in Section 2. In what follows, this summation will be understood in any formula of 
the emissivity. Contraction of the correlation function ( p^ with the neutrino-pair tensor KpR^ — K^gp^, as indicated 
in Eq. (p5|), yields: 



(K^K, - K^g^,) ^^-^ (K) = -K' 



1 



Im 



r<2 7^2 



47re2 exp (^) - 1 

T^AVT^VA 



+ Cl [ni + 2nt) 
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Taking the imaginary part of the right-hand side of this equation, after a lengthy (although straightforward) calcula- 
tion, we obtain the following formula : 



{K^K, - K^g^,) l.^'^ {K) = K' 
llm iTi\ 



47re2 exp (y) 



llm TTtl 



-Cl(|Im7r,| + 2|Im^t| 



+2Cl 



llm ttA 



(56) 



Here it is assumed that ImTr^ — ImTr^y — , and we incorporated the fact that the imaginary parts of retarded 
polarization functions are negative for uj > 0, so that — ImTT; t = \^'ccLT^i,t\- 

As the polarization functions Eqs. (|2^, |27|) and Eqs. ( ^ ^6|) are independent of the direction of the vector k, 
integration over angles can be done in Eq. (pa), and we obtain the vv emissivity from the plasma as follows 



487r5e2 



GlK^ 



dhj L 



dk- 



exp (^) - 1 
Im TTi I 



|Im TTf] 



RCTT;)^ + (ImTT;)'^ 



+Cl(|Im7rn +2|Im^t|) 
+^Cl—(ReTrAv) 



{K^ - Rent) + {limrtY 



|Im TTt\ 



K 



(57) 



This formula describes the vv emissivity from the plasma with inclusion of collective effects due to the plasma 
polarization. As we included the in-medium photon contribution to the correlation function, Eq. ( |57| ) describes the 
energy losses due to plasmon decay and vv Bremsstrahlung from electrons. Due to collective effects, the contribution 
to these particular processes through the vector weak-current can not be separated from one another. This can be 
understood as follows. The vector transition current of the electron is responsible for production of neutrino pairs and 
on-shell photons in the medium as well. When a virtual photon participating in the vv Bremsstrahlung goes on shell, 
it has such a large life-time, that decay of the photon into neutrino pairs can be dually interpreted as free plasmon 
decay, or as vv Bremsstrahlung from electrons. This interference does not permit to separate the partial contributions. 
However, such an interference can be neglected in some limiting cases, when one of the processes strongly dominates. 



VI. LIMITING CASES AND NUMERICAL TESTS 



A. Moderate temperatures 



At moderate temperatures T 



plasmon decay into neutrino-pairs strongly dominates the neutrino energy 



losses from the degenerate plasma of a liquid crust ]16| ] . In this case, the imaginary parts of the polarization functions 
result only in a minor widening of the plasmon spectra due to plasmon capture and creation in electron collisions 
with nuclei. By taking in Eq. ( |57| ) the imaginary part of polarization functions equal to zero, we obtain: 



where 



Q = Qi 



duj UJ 



exp 1 7^ 



Q 



t — 1 



Jo 



(58) 



(59) 



Qt 



G^p 



doj u 
exp (^) 



dk k^ {C^K^ + C\k'^TT\y) 5 {K^ - TTt) . 



(60) 
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As it follows from the poles of Eq. (49), longitudinal and transverse eigen-photon modes satisfy the dispersion relations 

cj2-fc2_Re7r;,t(cj,fc) = 0, (61) 



therefore Eqs. (59) and (|60|) are easily interpreted as the emissivity due to decay of on-shcU longitudinal (l) and 
transverse (i) plasmons, respectively, into neutrino pairs. We will now show that, in fact, Eqs. (^ ) and ( |60| ) give the 
known results for the neutrino pair emissivity from longitudinal and transverse plasmons, respectively. 
By expanding — Re ttj around ut one can approximate 



with 



uj^-k^- Re 7Ti,t {uj, k)~ — {uj^ - Luft {k)) 



(62) 



diTi^ti^, k) 



In this way 



Re Til-. 



2Lui,t{k) 



5{uj ~ uJLt{k)) , 



(63) 



(64) 



where LUi^t{k) is the solution (for > 0) of the dispersion equation (^l]) for longitudinal or transverse photons. 
Substitution into Eqs. (|59|) and (|60|) results in 



967r4e2 



duj 



1 



— 



dio I dkk 



exp ) - 1 7o 
' 2 Z,{k) 



exp ^ - 1 



dkk^ Zi{k)K^S{uj~uji{k)) 



C^K'^ + Clk^ReTTAvi^,k) 



5{uj - wt(fc)) 



(65) 



The integral over lo is trivial. As for the integral over k, one has to remember that the condition > implies that, 
for longitudinal modes, k can not be larger than the limiting value fcmax [0- Finally, we obtain 



Qi — 

Qt — — 



r<2 n'i 



G'^p 
48^2 



dk 



k^Ziik) 



{uf{k)-e)\ 



dk- 



exp (^) 
k^Ztik) 



exp 



(^) 



- 1 



(^k)-k'^f + C^k^ Re TT^y (cot {k) , k) 



(66) 



These formulae coincide with the known expressions pO[ | for neutrino pair emission from plasmons. 



B. High-temperature limit 

We now consider a degenerate electron gas in the case of high temperatures T 3> Wpe- Then the emissivity of 
Bremsstrahlung, which increases as T^, strongly overcomes the emissivity due to plasmon decay, which only goes 
as (see , e.g. [0). At such temperatures, the dominant contribution comes from ^ Keni^t ~ ^pe- Thus, 
neglecting Re7r;_t and Im7r;_f in the denominators of Eq. (p7|), and considering also (ReTr^y)^ ^ K^, we obtain 

Gi(ci + ci) dujuj r , , 

QBr= / dkk'K'{\lmn,\+2\lm7r,\). (67) 



48^5e2 exp(f)-l7o 

Performing this integral, with Imwi^t defined by Eqs. (|3^, ^6|) , we arrive to the known result ||ll| for the 
Bremsstrahlung emissivity: 



47r 



= 189^^ + e^A^.Z^T^ 



X 







^-rf.^(.)(^_2) r,,!i(i^l,C±-, (68) 
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The latter integral can be done analytically (see also H) : 



f\ u(l~u^) r + xu 2 1 / 2a;2 ^ \ 



(69) 



and we finally obtain 



0* = id + ca .««,z'r» [" ^ (i + JfL ,„ . ) . ,70) 



C. Low-temperature limit: collective effects 



In the case of low temperatures T <C ujpe, the occupation numbers of eigen-photon modes in the medium are expo- 
nentially suppressed, and the typical energy of a neutrino-antineutrino pair is of the order of the medium temperature. 
Under this situation we can neglect ^ with respect to Re 7r/^f ~ lo^^ in the denominators of Eq. ( ^7|) . 

Now, Eq. (57) takes the form 



(71) 



These terms describe the vector (Q^) and an axial (Q'^) weak-current contributions 



Q 



V — 



^2 /-i2 



exp (f ) - 1 



dkk^K^ 



VRe^ 



|lm TT/I 



if 4 



Re^TTt 



|Im TTt 



(72) 



In the latter equation we neglected terms which are proportional to ReTr^y/ReTTj ~ T/ ^ 1. 

As mentioned above, the imaginary parts of polarization tensors are due to electron collisions with nuclei. Thus, 
Eqs. (^, ^) should be understood as the neutrino-pair Bremsstrahlung from electrons. Substituting in Eq. ( [T^ ) 
the explicit form of the polarizations, and performing integrations over dk and duj we obtain, for a liquid crust, the 
following expression for the axial weak-current contribution 

Qi„ . '^ClCl.'N,Z^T' f" (l + J£L ,„) (74) 

567 7o \e{q)\q\ l~x^ J 

which coincides with the known result Eq. ( |70| ) for the axial term in the ultrarelativistic case, while the contribution 
of the vector weak current contains in the integrand the additional, small factors, / Re^ ni^t ~ T"^/^pe ^^'^ therefore 
decreases along with the temperature as T^°. Thus, at low temperatures, the contribution of the vector weak current 
to emissivity is suppressed due to collective effects, as discussed in 0. 

In Fig. 3, we show the result of a numerical calculation of the neutrino-pair emissivity from a liquid crust with 
a density p such that pY^, = 7.3 x 10^ gr cm~^ (Y^ is the electron fraction per baryon), consisting on nuclei with an 
atomic number Z = 30 embedded into a degenerate ultrarelativistic electron gas, which is in good agrement with 
the total neutrino emissivity, as given by the sum of the plasmon decay contribution, calculated in [|6|, and the 
"standard" Bremsstrahlung contribution calculated in ||ll| . One readily observes from this figure, that our calculation 
coincides with the standard Bremsstrahlung emissivity, when this is the dominant process (at high temperatures), 
and with plasmon decay at lower temperatures. 

Collective effects are dramatic at low temperatures T <^ ujpe, when plasmon decay is negligible, which are typical 
for the crystallin crust. However, the above criterion T > 172 for crystallization of the crust corresponds to the case 
where an infinite time scale is allowed to attain thermal equilibrium. In actual stellar evolution, only a finite time 
scale is allowed for temperature variation. Therefore, the ionic system is likely in an amorphous-crystal state for 
172 < r < 210 Thus, in practice, the liquid state of the crust persists up to F < 210. At such temperatures, the 
plasmon decay is exponentially suppressed, and Bremsstrahlung becomes again dominating. In this case, collective 
effects dramatically suppress the vector weak-current contribution to neutrino-pair Bremsstrahlung from electrons in 
the liquid crust. 
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Unfortunately, the structure factor is not tabulated for the liquid ionic system in this domain of T. For this reason, 
in order to make these effects more apparent, we considered a liquid crust with the same value of pKe, consisting on 
nuclei with Z = 10. Such a crust remains liquid for temperatures much smaller than the plasma frequency. The 
result of our calculation is plotted in Fig. 4, where we show separately the and pieces, in comparison with 
plasmon decay. The contribution to the emissivity shows the standard behavior ~ T^, while the piece drops 
much faster with temperature, as r^°. As we discussed above, such a behavior is a consequence of collective effects, 
which result on the screening of the vector weak-current of electrons in the plasma. 

In Fig. 5, we have plotted the ratio R of the total emissivity, as calculated from our formulae, to the sum of the 
standard Bremsstrahlung plus plasmon decay (taken from and |l^). As it readily seen from this figure, this ratio 
abruptly decreases at low temperatures, and tends to the expected value H. 



QO Cy + C\+2 {Cy 



VII. CONCLUSIONS 



We have shown a method of calculation for the neutrino-pair energy losses from the crust of a neutron star, which 
has the advantage that collective effects are automatically included, and do not need a separate consideration. The 
formula we obtained for the neutrino-pair emissivity, Eq. (^7|), includes the contribution of plasmon decay into 
neutrino pairs, as well as neutrino-pair emission due to electron collisions with nuclei and phonons. Some limiting 
cases, considered in Sect. 6, demonstrate that this formula describes the known results for energy losses due to 
plasmon decay and neutrino-pair Bremsstrahlung of electrons, when one of these processes dominates, at moderate or 
high temperatures. Our formula takes into account collective effects, which in the case of low temperatures manifest 
in a suppression of the vector weak-current contribution to the emissivity, in accordance with the results in |^ . Some 
numerical tests for moderate and high temperatures show a good agreement of our calculated total emissivity, with 
that obtained by summation of the separate contributions of plasmon decay and electron Bremsstrahlung obtained by 
different authors. In general, Eq. ( |57|) c an also be used for a non-degenerate plasma, if the corresponding polarization 
functions are inserted. Thus, Eq. (pTj) is able to incorporate neutrino-pair energy losses from the electron plasma 
caused by processes other than plasmon decay and neutrino-pair Bremsstrahlung. This can be done by improvement 
of the imaginary parts of the retarded polarization functions, with inclusion of different physical processes leading to 
damping of the plasmon in the medium under consideration. 
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FIG. 1. Diagrams contributing to the weak correlation function, when collective effects arc included. Shadow rectangles on 
electron loops represent intermediate electron interactions with rmclei (in the liquid crust) or with phonons (in the crystall). 
The thick dashed line on the second diagram stands for an in- medium photon (see next figure). 




FIG. 2. Series of diagrams defining the in-medium photon propagator (shown as a thick dashed line). Thin dashed lines 
correspond to the in-vacuum photon propagator. 
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T/w 

FIG. 3. Total neutrino-pair emissivity, calculated from our Eq. pA (solid line), in comparison with "standard" (i.e., with- 
out collective effects) plasmon decay (dashed line) and Bremsstrahlung (dashed-dotted line). All curves are for a density 
corresponding to pYe = 7.3 x 10^ grcm~^ and an atomic number Z = 30 for ions. 
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pe 

FIG. 4. Separate vector (solid line) and axial (dashed line) contributions to the neutrino-pair emissivity. At low temperatures, 
the vector weak-current contribution is suppressed with respect to the axial term due to collective effects. At moderate 
temperatures, when plasmon decay dominates, the vector weak-current contribution coincides with the plasmon decay curve 
(dashed-dotted line). All curves correspond the same value of pYe as in the previous figure, and an atomic number Z = 10 for 
ions. 
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FIG. 5. The ratio of the total neutrino-pair emissivity, as obtained from our formula, to the sum of standard plasmon decay 
and Bremsstrahlung. Plasma conditions are the same as in the previous figure. 
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